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i— i ABSTRACT 

Context. CO observations have been so far the best way to trace molecular gas in external galaxies, but in low metallicity environments 
the gas mass deduced could be largely underestimated due to enhanced photodissociation of the CO molecule. Large envelopes of FL 
j_( . could therefore be missed by CO observations. 

Aims. At present, the kinematic information of CO data cubes are used to estimate virial masses and trace the total mass of the 
molecular clouds. Millimeter dust emission can also be used as a dense gas tracer and could unveil H 2 envelopes lacking CO. These 
different tracers must be compared in different environments. 

Methods. This study compares virial masses to masses deduced from millimeter emission, in two GMC samples: the local molecular 
clouds in our Galaxy (10 4 -10 5 MO), and their equivalents in the Small Magellanic Cloud (SMC), one of the nearest low metallicity 
dwarf galaxy. 

Results. In our Galaxy, mass estimates deduced from millimeter (FIRAS) emission are consistent with masses deduced from gamma 
■ ray analysis (Grenier et al. 2005) and therefore trace the total mass of the clouds. Virial masses are systematically larger (twice on 

average) than mass estimates from millimeter dust emission. This difference decreases toward high masses and has already been 
reported in previous studies. This is not the case for SMC giant molecular clouds: molecular cloud masses deduced from SIMBA 
millimeter observations are systematically higher (twice on average for conservative values of the dust to gas ratio and dust emissivity) 
than the virial masses from SEST CO observations. The observed excess can not be accounted for by any plausible change of dust 
properties. Taking a general form for the virial theorem, we show that a magnetic field strength of ~ 15yuG in SMC clouds could 
provide additional support to the clouds and explain the difference observed. 

Conclusions. Masses of SMC molecular clouds have so far been underestimated. Magnetic pressure may contribute significantly to 
their support. 
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1 . Introduction structures dAndre et alj|2004l) . Most of the mass may then lie in 

the cloud envelopes where the gas is molecular due to H? self- 

Molecular clouds are of particular importance as the sites of star shielding but co is ph otodissociated dIsrae¥2005clLeWux et all 

formation in a galaxy. Most of this gas is located in large com- Q 994) Millime ter dust emission could bring to light the total 

plexes known as giant molecular clouds (GMCs), but the ex- mass f sucn clouds 
act amount of gas contained remains unknown since H2 is quite 

impossible to observe directly in cold interstellar regions where The advent of sensitive bolometer arrays on sub- 
molecules form and survive millimeter/millimeter telescopes makes long wavelength 
Carbon monoxyde (CO) has played a major role in GMC dust emission an alternative to CO to compute molecular 
studies and is still the most widely used tracer of molecular masses - The 8 ood correlation with CO emission in CO 
gas. However, it could be present only in the densest parts bright galaxiesje.g. NGC891 | Guelin etalj dl993|), NGC6946 
of molecular clouds, leaving H 2 envelopes unseen. This effect l Alton et aL l 'I™)) asserts lts a ssociation with molecular gas. 
has been clearl y observed in Molecular Clouds of the Solar More generally existing observations show that sub-millimeter 
Neighbourhood dGrenier et al.ll2005l) and tends to underestimate dust emission from galaxies traces the total hydrogen column 
molecular cloud masses estimated from CO observations. This denslt y (Thomas etalj | 2002t i Menerink et alj | 2005|) . It comple- 
bias is most important in low metallicity environments where ments the Hl and CO observations to provide a complete view 
the reduced shielding causes CO to be confined to tiny dense of cold interstellar matter. 

We propose to use dust continuum emission at millimeter 

Send offprint requests to: C. Bot wavelengths to compute molecular clouds masses in two sam- 
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pies of GMCs: molecular clouds in the Solar neighbourhood, and 
counterparts in the Small Magellanic Cloud (S MC). The SMC 
is a ne arby irregular dwarf galaxy (D« 61kpc: IL"anev &_S tobie 
1994)) known to have a low metallicity (Z^Zn/10. iDufou" 
1984)) and giant star forming regions. While the amount of 
molecular material observed via CO emission is low, numer- 
ous young stars are observed. Furthermore, galaxy counts be- 
hind the SMC imply extinctions much large r than expecte d from 
the column density of atomic hydrogen (Lequ euxlll994l) . Thus 
the SMC is a good candidate to detect extended envelopes of H2 
gas through dust observations. 

From the comparison of GMCs CO luminosity and virial 
mass, the SMC Xco = N{H2)/Ico conve rsion factor is esti- 
mated to be 10 tim es the Galactic value (Ru bio et all [1993a; 
Mizuno et al. 2001). This virial calibration of Xco should take 
into account the cloud mass in the CO-poor H 2 gas, b ut dust ob- 
servations question this view. iLisenfeld et alj ((2002) inferred a 
Xco conversion factor 25 to 30 times higher than the Galactic 
value by comparing sub-millimeter dust and CO emission from 
the low metallicity dwarf galaxy NGC1569. iRubio et al.1 f2004) 
reported the detection of the millimeter dust emission from 
SMCB1#1 with SIMBA/SEST. They converted the flux into a 
cloud mass estimate that is found to be one order of magni- 
tude larger than its CO virial mass. lLeroy et al.l (12006) analysed 
Spitzer far infrared maps in the SMC, in conjunction with Hi and 
CO observations. Comparing the 160/im map to NANTEN CO 
observations they estimate Xco to be 60 times the G alactic value. 
This p resent study follows and continues the work of Rubi o et all 
(2004) by analysing a larger number of molecular clouds in the 
SMC and by analysing these clouds together with a reference 
sample that we build in our Galaxy. As such it is the first compre- 
hensive comparison of two giant molecular cloud mass estimates 
at different metallicities: gas masses from dust millimeter emis- 
sion and virial masses. Are the dust and virial mass estimates 
in agreement? If they are discrepant, how can the difference be 
explained? 

After an overview of the data used in this study (Sect. [2]), 
molecular gas masses from dust millimeter emission are com- 
puted for a sample of Local GMCs and for detected SMC GMCs 
(Sect. 13. U . A comparison to virial mass estimates is performed 
in Sect. l3.2l and the differences are interpreted in Sect. [4] 



Table 1. 

SMC. 



Observed characteristics of the detected GMCs in the 



2. The data 

This study follows the detection of several molecular clouds at 
1.2 mm in the SMC with the SIMBA bolometer. We present 
these new observations in details, along with previously pub- 
lished data. 



2.1. SEST observations in the SMC 
2.1 .1 . CO observations 

The CO observations of the SMC clouds were observed 
as part of the ESO/S EST Key Prog r amme on CO in th e 
Mage llanic Clouds (llsrael et all 11993b IRubio et alJ ll993allR 
1996; Isra el et aT]|2003l) . The motivation for undertaking these 
observations in the transitions J=(l-0) and J=(2-l) of the CO 
molecule, and the reduction treatments are described in these 
papers. We just recall here the spatial resolution: the beam has a 
FWHM of 43"and 22"at the CO(1-0) and the CO(2-l) frequen- 
cies respectively. Since CO(2-l) line data has better resolution, 
gas masses computed for the CO(2-l) emitting regions will be 
favored when CO(2-l) observations were performed. In the fol- 



name 


Area 


S 1.2m 


a 
>: _ 




pc~ 


mjy 




LIRS49 


357 


181 


±55 


LIRS36 


385 


184 


±48 


SMCB1#1 


108 


47 + 


7 


SMCB1#2 


260 


119 


±51 


SMCB1#3 


38 


5.8 + 0.7 


Hodge 15 


103 


40.4 


±2.5 


SMCB2 S 


695 


398 


± 158 


SMCB2 N 


314 


132 


±46 


Nil 


87 


34.5 


+ 5 



corrected from the contribution of free-free emission and CO(2-l) 



line 



lowing, we will use the spatial (radius) and kinematic (velocity 
dispersion) characteristics of the clouds that were computed in 
the litterature. 



2.1.2. SIMBA millimeter continuum observations 

The SIMBA bolometer was mounted at the Swedish ESO 
Submillimeter Telescope (SEST). It operated at a central fre- 
quency of 250 GHz (1.2 mm) with a FWHM bandwidth of 90 
GHz, and an angular resolution of 24" on the sky. The indi- 
vidual maps were produced using the fast scanning mode. The 
elimination of the correlated sky noise, adding of the individual 
maps, and the photometry were done using the MOPSI pack- 
agqj. The T zen \th was determined performing skydips on average 
2 hours during the day, and every 3 hours during the night. The 
flux determination was done for each individual observing run 
by taking maps of Uranus. The data reduction was performed as 
described in the SIMBA User manuaQ. The flux measurements 
in the final coadded map has an uncertainty of 15%. 

The millimeter dust emission in the SMC is faint and 
the regions unambigously detected are giant molecular clouds 
(GMCs) (see FigQ]). Mo st of them have been o bs erved in CO an d 
are named according to Rubi o et ail dl993albl) or lHenizel dl956). 
SIMBA fluxes for each region are displayed in TableQ] They are 
measured for a SIMBA beam convolved to 27" and within the 
25 mJy/beam contour, after subtraction of the flux measured in 
an empty zone. 

Detected sources: Detected GMCs are often associated with 
well-known H 11 regions and you ng clusters an d all but LIRS49 
correspond to absorption patches (Hod gell974l) . We review here 
some of their properties: 

- LIRS49 and LIRS36 are tw o of the most prominent SMC 
sources in the far-infrared dSchwering & Israeli 1 19891) and 
a large variety of molecular specie s have been observe d 
in these directions dChin et alJll998t iHeikkila et alj fl998). 
LIRS49 may be a superposition of two H n regions belong- 
ing to different kinem atical complexes along the line of sight 
(Rosa doet al.lll994|) or it could be as sociated with a bubble 
or a loop (iMeyssonn ier & Azz opardi 19931). 

- SMCB1#1 and Hodgel5 are two quiescent molecular clouds 
(not associated with Ha). The SIMBA millimeter observa- 



1 MOPSI is a data reduction software package for IR and radio data 
developped by R. Zylka, IRAM, Grenoble, France 

2 http://puppis.ls.eso.org/staff/simba/manual/simba/ 
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Fig. 1. SIMBA mosaic of the south of the SMC (left panel). Names of the star forming regions studied in this paper are labeled 
closeby. On the right panel, CO contours of the studied regions are displayed on the same grid (with plain lines when CO(2-l) is 
available, with mixed lines for CO(l-O) data). 



tions of SMCB1#1 were used by iRubio et al] d2004l) to de- 
duce a cloud mass. 

- SMCBW2 and SMC B1#3 are parts of the SMC-B1 region of 
Rubioetal. (1993b]) and are associated with large H n re- 
gions. 

- The SMCB 2 region contains several molecular clouds ob- 
served by Ru bio et all {1993a). Since these clouds are diffi- 
cult to separate in the continuum dust emission which has 
no velocity information, they were grouped into two com- 
pl exes: SMCB2 North (corresponds to SMCB2#3 and #4 
of IRubio et all dl993al) ) and SMCB2 South (which corre- 
sponds to SMCB2# 1, #2 and #6). The SMCB2#5 region of 
IRubio et al]dl993al) is not detected at 1.2mm with SIMBA. 

- Millimeter e mission is also detected toward the N30 and Nil 
dHenizd [T956) H ii regions, but these regions have not yet 
been mapped in CO, so no mass comparison can be done for 
these objects. For N30, the SIMBA flux can be accounted for 
entirely by free-free emission. 

- The source identified by the question mark in Fig. Q] (be- 
tween SMCB2 North and South) does not correspond to any 
presently known source. It is clearly detected at the 6cr level 
but there is no associated infrared or radio component. It 
should be noted however that it lies on an extended compo- 
nent seen at 70/im with Spitzer observations and it is close 
to the DEM S43 H n region. Further observations are needed 
to understand the nature of this source. 



The observed regions correspond to different evolutionary 
stages of molecular clouds - from quiescent clouds to evolved 
massive star forming regions-, enabling us to compare CO and 
dust mass estimates in different physical conditions. 



Origin of millimeter continuum emission: The emission mea- 
sured at 1.2 mm consists of several components: thermal dust 
emission, free-free radiation and 12 CO(2-l) line emission. Since 
our interest is only in cold dust emission, we estimate the other 
contributions that have to be removed. 



CO line emission: The I2 CO(2-l) line contribution to the 
SIMBA surface brightness can be estimated as follows: 

2/fcvV 3 

Fline - — IcO(2-\), (1) 

Avbol 

where Ico(i-\) is the intensity of the I2 CO(2-l) line in K.km.s" 1 
, Av ho/ is the bandwidth of the SIMBA bolometer (90 GHz), 
Q, is the region of integration, k is the Bol tzmann constant an d 
v = 230 GHz. For regio ns where CO(2-l) dRubio et al.ll 1993bllaL 
ll996Hlsrael et al.l20 03) has been observed, the total flux density 
due to line contributions was evaluated (see Table [TJ and was 
found to be negligible. 

free-free emission: The contribution of the continuum emis- 
sion due to thermal electrons can be estimated with radio contin- 
uum data, assuming pure free-free emission at high frequencies 
and extrapolating this frequency dependence to 1.2 mm (250 
GHz) B We also assume that radio emission is optically thin 
at long frequencies. W e used the ATCA continuum catalog of 
the SMC radio sources dPavne et al.ll2004l) at 1.42, 2.37, 4.8 and 
8.64 GHz. In both SIMBA and ATCA observations, extended 
emission on scales larger than a given threshold is filtered out. 
For SIMBA, extended emission on scales larger than ~ 2' is fil- 
tered out when subtracting the sky emission. The interferometric 
ATCA observations are not corrected for missing short spacing 
and do not measure emission on spatial scales larger than about 
half the size of the primary lobe. At 8.64 GHz, the two angular 
scales are comparable. Radio emission comes both from free- 
free and synchrotron emission but at these high frequencies it is 
dominated by free-free. We estimate the free-free contribution 
by extrapolating the 8.64 GHz radio flux when available, since 
it is the less synchrotron biased frequency of this catalog and 
also the resolution is the closest to the SIMBA one. When no 
8.64 GHz flux is relevant, we take the longest frequency avail- 
able and follow the same procedure: for LIRS49, the free-free 
contribution is extrapolated from the 2.37 GHz flux; for SMCB2 

3 For electron temperatures typical of Hn regions (~ 10 4 K ), the free- 
free e mission scales as (0.95 - 0.16 In v) where v is in GHz (Reynolds 
ri992h 
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North and South region it is obtained from the 4.8 GHz fluxes. 
It should be kept in mind in these cases that the free-free con- 
tribution could be overestimated, a bias that would reinforce the 
result (see section|4|. 

dust emission associated with neutral atomic gas Hi. In the 
SMC, the resolution of Hi observations does not allow an esti- 
mate of the contribution of the SIMBA flux that is associated 
with neutral gas. However, most of this emission is spatially ex- 
tended and filtered out in the subtraction of the sky emission. 

For each source, dust emission is the major component of the 
1 .2mm flux. In the following, we shall refer to the SIMBA flux 
tabulated in Table [T] as the flux measured at 1.2 mm subtracted 
from the gas contribution estimates. 

2.2. FIRAS and CO data in our Galaxy 

To complement the analysis of SMC clou ds, we use F I RAS spec- 
tra and the CO (1-0) sky survey from IDame et all (1200 ll) for 
a sample of Local GMCs defined in rectangular (/, b) boxes^: 
Orion, Taurus, Ophiuchus, Cepheus, Chamaeleon, Lupus and 
Perseus. The 7" FIRAS beam for these clouds at a few hundreds 
of parsecs and the SIMBA beam at the distance of the SMC cor- 
respond to similar linear spatial resolutions. 

2.2.1. Millimeter fluxes 

The 1 .2mm fluxes were derived by fitting the FIRAS spectra at 
A > 500/im with a black body multiplied by a power law emis- 
sivity. The temperatures in Table [2] were derived from a single 
temperature fit of the full FIRAS spectra down to 150/jm and 
an emissivity index of 2. As for the SIMBA observations, the 
1.2mm emission observed toward the Galactic molecular clouds 
originate from different components other than cold dust emis- 
sion. We try to estimate their contribution to the 1.2mm flux. 

The CO(2-l) line emission contribution is difficult to es- 
timate without CO(2-l) data on hand. However, since bright- 
nes s temperatures are the same in CO J=2-l and J=l-0 lines, 
the IDame et all (1200 ll) survey can be used to estimate this con- 
tribution. Usin g a Ico<.i=2 7 \\l Ic oli=\-Q) - 0.6 for the Solar 
Neighborhood ( Hasegawal 1 1 9961) . this contribution is found to 
be negligeable. 

Free-free emissio n: extinction corrected-Ha emission 
dDickinson et al.l2003l) is used to estimate the free-free emission 
contribution to the FIRAS fluxes. This is found to be very low 
(< 3%) and is therefore negligeable. 

Dust emission associated with neutral atomic gas. Unlike 
the SIMBA observations, the FIRAS data measures the to- 
tal Galactic emission. Part of this emission comes from 
dust associated with the diffuse medium. We used the 
Leide n/Argentine/Bonn (LAB) Hi survey of Kalberl a~et al.l 
(2005) to estimate the proportion of millimeter emission asso- 
ciated with neutral gas. At high galactic latitude, the dust mil- 
limeter emission is w ell correlated with Hi column densities 
dBoulanger et alj[l996l) . The dust emission per H atom is com- 
bined with the Hi column densities measured in the cloud areas 
to compute and subtract the 1.2mm flux associated with neu- 
tral gas. Their contribution amounts to about 40% of the total 
FIRAS flux at 1.2mm. The 1.2 mm fluxes obtained are listed 
in Table [2] for each cloud. To estimate the uncertainties associ- 
ated with the dust emission removal, we defined for each cloud 
a reference area of similar size on the symetric position within 



4 The box were defined on the millimeter emission map 



the Galaxy (same longitude, opposite latitude) and subtract the 
flux observed in the reference regions from the one observed in 
the molecular clouds. The differences observed are similar to the 
uncertainties on the data quoted in Table|2] 

2.2.2. CO emission 

We use the CO(1-0) survey of IDame et all (120011) of the Galaxy 
to compute the size and velocity dispersion of each cloud over 
the same areas as the millimeter emission observations. To be 
consistent with CO information in the SMC, we chose to com- 
pute the size and the velocity dispersion of each cloud in a simi- 
lar way (e.g no correction for velocity gradients). For each cloud, 
we fit each spectrum above a certain level and sum the gaus- 
sian fits to deduce a total cloud spectrum. The velocity disper- 
sion of the cloud is determined by fitting a gaussian to the re- 
sulting spectrum. This method is close to what was used for 
the SMC except that some SMC clouds were divided into sub- 
structures when independent clouds were separated in spatial- 
velocity space. We checked that this difference in method for 
some clouds does not induce a significant bias in the cloud pa- 
rameters. The radius and velocity dispersions are listed in Table 
[2] for each cloud. Because there is insufficient spatial resolution 
to correct for velocity gradients in the SMC clouds, we chose not 
to correct for them in the Galactic sample either. We are careful 
to apply the same method to deduce virial masses in both data 
sets. 

3. Estimates of molecular cloud masses 

3.1. Molecular gas mass from millimeter continuum emission 

Dust millimeter continuum emission is associated with molecu- 
lar gas. In this section, we compute molecular cloud masses from 
the dust emission at 1.2mm. 

The far-infrared/millimeter emission from dust can be ex- 
pressed as: 

h = N H e H (A)B A (T dusl ) (2) 

where Nh is the hydrogen column density of the cloud, en{A) 
is the emissivity of dust per hydrogen atom at the wavelength A 
and T c iust is the dust temperature. The emissivity per hydrogen 
atom can be expressed in terms of the absorption coefficient per 
unit dust mass k(A): 

e H (A) = K(A)xdHm H (3) 

where xj is the dust-to-gas mass ratio, and ymiu is the gas weight 
per H, taking into account the contribution of He. This emissiv- 
ity may be determined empirically from the correlation between 
dust and gas obs ervations, independently of any dust model (e.g 
iBoulanger et al.l (119961) for the diffuse interstellar medium in the 
Solar Neighborhood). 

While dust temperatures for the Galactic clouds have been 
deduced from the FIRAS spectra fitting , this information is not 
available for the clouds in the SMC. Since we wish to apply the 
same method for both samples, we will fix the dust tempera- 
ture for both samples. Within each cloud, dust is present over 
a range of temperatures. The fit of FIRAS spectra with a sin- 
gle temperature gives higher weight to the warmer dust present 
at the cloud surface or near embedded stars. But in a molecular 
cloud, the bulk of the mass is shielded from UV radiation. The 
dust temperature in the shielded parts of the clouds is observed 
to be smaller than the values in Table [2] For example, based on 
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Table 2. Observed characteristics of the Galactic GMCs studied. 



name 




distance 






R 


AV 






pc 


mjy 


K 


pc 


km.s 


Orion a 


204.6, 216.4/ -23.9, -13. 6 


450 


11000 ± 3000 


19.1 


28.9 


8.0 


Orion b 


202.1, 205.9/ -14.4, -11.1 


450 


1700 ± 300 


19.1 


11.6 


3.6 


Taurus 


170.1, 188.9 / -19.9 , -10.1 


140 


14000 ± 4000 


16.2 


13.5 


3.8 


Ophiucus 


347.1,359.9/ 11.6,24.9 


130 


15000 ± 4000 


19.6 


11.5 


3.8 


Cepheus 


102.1, 116.9/ 12.1 ,20.9 


300 


9700 ± 4000 


16.1 


22.5 


6.8 


Chamaeleon 


292.1, 305.9/ -19.9, -10.6 


160 


6100 ± 3000 


14.8 


6.2 


4.6 


Lupus 


330.1,344.9/ 10.1 , 19.9 


140 


7000 ± 4000 


17.4 


5.2 


5.1 


Perseus 


156.6, 161.9/ -22.9, -15.1 


320 


4700 ± 1500 


17.1 


14.6 


6.9 



" corrected from the contribution of dust associated with H i 



submillimeter observations, Stepniketal. (2003) report a dust 
temperature of 12 K for a Taurus filament. To conservatively es- 
timate cloud masses, we adopt a dust temperature of 15 K (the 
minimum value of the FIRAS temperatures listed in Table |2]i 
to derive Galactic and SMC cloud masses from their millimeter 
dust fluxes. 

We assume that the dust grain opacity per unit mass is the 
same for the SMC and Milky Way molecular clouds. Therefore, 
the emissivity per hydrogen atom scales as the dust-to-gas mass 
ratio. Knowing the SMC dust-to-gas mass ratio and the Galactic 
dust emissivity per hydrogen atom in the SMC, the millime- 
ter emission can therefore be translated into a gas amount. The 
Galactic emissivity per hydrogen atom is adopted for dust asso- 
ciated with molecular gas (see Appendix lAl: e r ^ (1.2mm, H2) = 
(1.2 + 0.2) x \Q- 2(, ar l .cm 2 . For the SMC, this emissivity 
will be reduced by one sixth to account for the difference in 
dust-to-gas ratio with respect to the Solar Neighbourhood (see 
Appendix|B]). The values of dust emissivity per hydrogen atom 
we use, correspond to absorption coefficient per unit mass of 
k = 0.72cm 2 .g~ 1 , 3 times higher than the standard value from 
iDraine &~ Lee (1984). Furthermore, the du st-to-gas ratio i n the 
SMC could be as much as 5 times smaller dBot et alj|2004f) . The 
mass estimates we deduce from millimeter emission in the SMC 
are therefore upper limits. These choices were made to empha- 
size the robustness of the results of this study, but the range of 
parameters available for the dust emissivity and dust-to-gas ratio 
is assessed throughout the paper (see Sec.|4]and Fig. [2]). 

We can thus translate the measured millimeter fluxes of the 
clouds into their corresponding hydrogen column densities and 
molecular cloud masses. The results for the Solar Neighborhood 
and the SMC clouds are presented in tablesOandlHrespectively. 
Note that we used an emissivity per unit dust mass and a dust- 
to-gas ratio tha t are 1.8 and 1.7 tim es higher, respectively, than 
the one used bv lRubio etail d2004l) . 

3.2. Comparison to virial masses 

In both data sets, the kinematic and spatial information provided 
by the CO data cube enables the determination of virial masses. 
For this determinati on, we use the simplifie d formulation for a 
1/r density law from lMacLaren et akl (11988): 

M vir (M ) = l90AV 2 (km.s- [ )R(pc) = (4) 

where AV is the observed velocity dispersion from the CO 
linewidth (FWHM) and R is the radius of the cloud. The virial 
masses for both cloud samples are computed using the radius 



Table 3. Different mass estimates of the Galactic GMCs studied. 



name 


10 4 M o 


AT"" 

H 

10 4 M o 


/ = 


T 

M H 


Oriotf 


14 + 3 


27 + 2 


0.5 


±0.2 


Taurus 


1.3 + 0.4 


2.6 + 0.1 


0.5 


±0.2 


Ophiucus 


1.3 + 0.3 


2.3 + 0.1 


0.6 


±0.2 


Cepheus 


4.3 + 1.8 


14.2 + 0.8 


0.3 


±0.2 


Chamaeleon 


0.8 ± 0.4 


1.7 + 0.1 


0.5 


±0.3 


Lupus 


0.7 ± 0.4 


1.8 + 0.1 


0.4 


±0.3 


Perseus 


2.4 ± 0.7 


9.3 + 0.4 


0.3 


±0.1 



" the mass are obtained by adding the masses of Orion a and Orion b 

Table 4. Different mass estimates of the SMC GMCs studied 
with the associated uncertainties. Notes: (1) computed from R 
and AV using CO (1-0) data; (2) computed from R and AV using 
CO (2-1) data 



name 


M""" 
H 

10 4 Mo 


AT"' 
10 4 M G 


1 < 


LIRS49 


19 ±6 


8.9 ± 0.5 (2) 


2.2 ± 0.7 


LIRS36 


19 ±5 


14.0 ± 0.5 <2) 


1.4 ±0.4 


SMCB1#1 


5.1 ±0.8 


2.5 ± 0.2 <2) 


2.0 ±0.3 


SMCB1#2 


13 ±6 


7.3 ± 1.4W 


1.8 ±0.8 


Hodge 15 


4.3 ±0.3 


4.9 ± 0.4 (2) 


0.9 ±0.1 


SMCB2 S 


43 ± 17 


14.3±2.3 (1) 


3.1 ± 1.2 


SMCB2 N 


14 ±5 


9.8 ± 1.9 (1) 


1.5 ±0.5 



and velocity dispersion of the clouds. For the SMC, these pa- 
rameters were taken from the litterature (c.f. Sec. 12. 1.1] ). For the 
Galactic clouds, the radius and velocity dispersion were com- 
puted as described in Sec. 12.2.21 in order to be consistent with 
the SMC values. The results for Galactic and SMC clouds are 
pre sented in tab les [3] and|4]respectively. 

iRosolowsk v & Lerovl d2006l) have shown that spatial resolu- 
tion and sensitivity differences in molecular-line data cubes can 
significantly affect the measured sizes and line widths that are 
used to deduce virial masses of giant molecular clouds. In the 
present study, the linear spatial resolution of the S EST CO(2- 
1) dat a in the SMC and of the CO(1-0) data from lDame et all 
(2001) in the solar neighbourhood are similar. Bias arising from 
low resolution should therefore affect the SMC and Galactic 
sample in the same way. This offset may be more problematic 
when no CO(2-l) data was available for some clouds in the SMC 
and CO(1-0) data with larger resolution (~ 13pc) was used in- 
stead. However, the results we obtain in these cases are not sig- 
nificantly different that the results obtained with higher resolu- 
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Fig. 2. Evolution of the ratio of the virial mass and the mass 
deduced from millimeter emission with respect to the mass de- 
duced from the millimeter dust emission in solar masses M .The 
Galactic clouds are represented by rhombs while th e SMC 
clouds are triangles. The law of ISolomon et al.l (Il987l) is also 
plotted (the solid line), representing the evolution of a larger 
sample of Galactic clouds. For the SMC clouds, the uncertain- 
ties on masses deduced from millimeter emission (regarding to 
the assumptions on the dust to gas ratio and the dust emissivity) 
are represented by a straigth line for each cloud. 



tion CO data so the uncertainties induced do not affect the results 
of this paper. 

Virial masses for the clouds of the Galactic sam- 
ple have been publish ed prev iou sly ( de Geus et al.l 
Ungerechts & Thaddeusl Il987t [Boulanger et alJ 
Grenier et al.l 1 1 989t> . However, the characteristics of 



1990; 



the 

clouds were computed for substructures and the masses of 
the substructures were then added. We have compared our 
virial mass estimates with those of the litterature for the same 
clouds. We observe minor differences that reflect the different 
definitions in the extent of the clouds and the difference in the 
method used (velocity gradients, substructure divisions). 



4. Results 

Tables [3] and |4] give mass estimates for the clouds in our Galaxy 
and in the SMC, computed from millimeter observations of dust 
emission and from virial estimates. These two estimates of the 
gas mass of a cloud are compared through the ratio of these 



A/„, 



masses, / = 

We observe a systematic difference between the two mass 
estimates: in our Galaxy, the molecular cloud masses deduced 
from millimeter emission are systematically lower than the 
virial masses obtained from CO observations, while in the SMC 
the mass deduced from millimeter emission are systematically 
larger. This difference is best seen in Figj2] 



4.1. In the Galaxy 

The comparison of the virial mass and the gas mass deduced 
from millimeter emission in the local Galactic clouds shows that 
the virial mass is always larger than the mass from millimeter 
emission, with a mean / factor of 0.5. 



iGrenier et al.l (|2005 ) analysed y ray emission in Giant molec- 
ular clouds of the Solar Neighborhood and found that a signifi- 
cant part of the mass of large molecular clouds is not traced by 
CO emission. We compared our mass estimates to their masses 
of "dark gas" (interpreted as H2 without CO) for the same re- 
gions. We find that the cloud masses deduced from millimeter 
dust emission are fully compatible with the total masses of the 
clouds (the "dark gas" plus molecular gas traced by CO emis- 
sion). However, for the regions chosen in this study, the fraction 
of "dark gas" with respect to the total molecular gas is low. 

ISolomon et alJ d 1987b compared, for a large sample of giant 
molecular clouds in our Galaxy, the virial mass of the clouds to 
their CO luminosity. They observed a good correlation between 
both quantities that can be transformed into a relation between 
the mass deduced from CO emission and the virial mass. Since 
in the studied regions the masses deduced from CO emission 
are close to the masses deduced from millimeter dust emission, 
this law can then be compared to our mass estimates for the 
Galactic clouds (see Fig. [2). For the Galactic clouds, the mil- 
limeter and CO masses are consistent for the Galactic gamma 
ray calibration of Xgp = 1.8 ■ 1 20 m ol.cm ^^K.km.s 1 ) 1 . For 
this factor, the Solomo n" et al.l (1 19871) virial mass-Leo relation 
show that the virial mass is larger than the cloud mass by a fac- 
tor decreasing toward high masses (see Fig. [2ji. This trend was 
also observed for clumps within GMCs by numerous studies 
in our Galaxy dBertoldi & McKed 119921; IWilliams et all 1 19951; 
iHever et al J [20011: ISimon et al.ll200ll) . lit the Galaxy, the / fac- 
tor is observed to be 1 only for the most massive Ga lactic GMCs 
with masses larger than 5 TO 5 M (iDame et al.ll986l) . For clumps 
within large gravitationally bound clouds the weight of the over- 
lying material in the cl oud could provide signi ficant binding 
pressure for the clumps (Bertoldi & McKee 1992). 



1998; 4.2. In the SMC 



The comparison of gas masses obtained from the millimeter dust 
emission with the virial masses shows a systematic bias in the 
SMC. Gas masses deduced from the millimeter emission are 1.8 
times higher on average than the virial masses. The factor / = 
■ ranges from 0.9 to 3. 1 (see Table|4|. The mass discrepancy 



reported for SMCB1#1 by iRubio etail (12004 is thus a general 
result that must be understood. 

For the computation of masses from the millimeter flux, we 
used a dust temperature of 15 K for all the clouds. Could a 
higher temperature in the molecular cloud of the SMC explain 
the difference observed in mass? For agreement between the 
virial mass and the mass from millimeter dust emission to agree, 
the dust temperature within the molecular gas would have to be 
23 K on average. But if the molecular clouds of the SMC were 
warm, they would emit in the far infrared. While for most of the 
regio ns surveyed here, the Spitzer 160^m emission dLeroy et al.l 
2006) is biased toward the Hn region emission present in the 
beam, this is not the case for the quiescent molecular cloud 
SMCB1#1. Comparing millimeter fluxes with Spitzer data in 
SMCB1#1 gives a dust temperature of 11, 14 or 17 K for spec- 
tral indices of 2, 1.5 and 1, respectively. This shows that at least 
in a quiescent cloud (and presumably in all molecular clouds) a 
higher dust temperature in the molecular clouds of the SMC is 
not compatible with the far-infrared data and can not explain the 
difference observed between the mass estimates. 

Furthermore, the cloud masses deduced from dust millime- 
ter emission depend on the dust-to-gas ratio and the emissivity 
of dust grains in the SMC molecular clouds. The values we used 
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for these parameters were chosen to limit the mass of the clouds 
and thus reduce the discrepancy with virial masses. In the diffuse 
ISM of the SMC, the dust emissivity per hydrogen atom is mea- 
sured to be 1/3 the emissivity of the Galactic diffuse medium 
dBot et all 2004). For the molecular clouds of the SMC, we used 
a larger emissivity to account both for grain coagulation effects 
(twice the diffuse medium value) and for our estimate of the 
SMC dust-to-gas ratio (twice larger than the canonical value for 
the metallicity of the SMC, c.f. Appendix iBli. Not only it is not 
clear if grain coagulation (grain emissivity enhancement by mass 
unit) occurs in the SMC, but we assume that heavy elements de- 
pletions are the same in the SMC and Galactic molecular clouds 
while they are observed to be smaller in the SMC diffuse ISM. 
In that sense, we consider the millimeter mass estimates in Table 
[4] to be lower limits on the molecular cloud masses. Lower de- 
pletions and thereby lower dust emissivity per hydrogen atom 
will enlarge the mass discrepancy. The range of masses obtained 
for each cloud, for dust emissivities ranging from 1/6 to 1/20 the 
Galactic molecular value, is shown in Fig. [2] Since it is dynami- 
cally difficult to account for a large difference between the cloud 
mass and its virial estimation, a first conclusion of this paper is 
that SMC dust observations imply a significant (up to a factor 
~ 5) increase in the dust-to-gas mass ratio from the diffuse ISM 
to molecular clouds. The dust which is destroyed by shocks in 
the diffuse ISM is rebuilt in molecular clouds. A second con- 
clusion is that even with conservative assumptions on the dust 
emissivity per hydrogen atom, the virial mass is found to be sys- 
tematically lower than the millimeter gas mass. 

5. Generalized virial theorem 

In the SMC, molecular cloud masses from millimeter dust emis- 
sion are systematically larger than the virial masses. We saw 
above that uncertainties on the mass estimates from dust emis- 
sion can not explain such differences. Observational estimates of 
the velocity dispersion are obtained through the CO line width. 
The lack of resolution and sensitivity of the observations cer- 
tainly plays a role on the virial mass computations. However, the 
observational bias on the Galactic and the SMC sample should 
be similar and we were careful to apply the same analysis, so we 
are confident that the virial mass discrepancy in the SMC does 
not a rise from o bservational uncertainties. 

iLerov et all d2006l) studied dust emission in the far-infrared 
in the SMC and showed that excesses with respect to Hi are 
observed to extend farther than CO emission observed with 
NANTEN. This difference is attributed to molecular hydrogen 
not traced by CO and amounts to 30% in spatial coverage. The 
CO and SIMBA observations we use for the SMC clouds cer- 
tainly underestimate the size of the giant molecular clouds due to 
sensitivity effects. As a result, the masses measured in this study 
are characteris tics of clumps w ithin the clouds smaller than the 
one from the ILerov et all d2006l) study. However, the size of the 
clouds that are detected with SIMBA and CO observations (see 
Fig- 03 in our study are about the same, so that the effect of any 
biases in the cloud radius should not be large enough to affect 
our results. 

The mass discrepancy observed between the SMC and 
Galactic clouds might indicate that the CO line width underes- 
timates the amplitude of the cloud turbulent mo tions. This in- 
terpretation was proposed by iRubio et"aT] d2004l) . In the SMC 
clouds, the CO emission is expected to come from dense clumps 
embedded in lower density gas where CO is photo-dissociated 
dLequeux et al.fl994l) . The bulk of the cloud mass is in the lower 
density gas but the CO line width only measures the velocity dis- 



persion of the clumps. The same statement applies locally to the 
less massive high Galactic latitude clouds where CO line widths 
are observed to be systematically lower than the Hi line width. 
Quantitatively, in the Ursa Major high latitude cloud, on scales 
of 5-10 pc, the Hi velocity disp ersion is observed to be twice 
that of CO dde Vries et alJll987l) . It is not clear to us whether 
this observed difference can be simply extrapolated to the more 
massive SMC clouds that are gravitationally bound. However, 
it has been theoretically suggested that in gravitationnal bound 
clouds the gas velocity dispersion shou l d increase as the den- 
sity decreases dFatuzzo & Adamsl [19931: lMcKedll999l) . In this 
section, we present an alternative interpretation where the mass 
difference results from the contribution of the magnetic field to 
the clouds support. 

In this section, we inspect the assumptions made for the 
virial mass calculations, searching for those which could cause 
a systematic underestimation of the mass. In particular, the pres- 
ence of a magnetic field in the molecular clouds could help to 
balance gravity and we will compute the necessary strength to 
explain the mass ratio / observed in the SMC. 

To apply the virial theorem, molecular clouds should be viri- 
alized. Observationally determined dynamical parameters en- 
able us to define the characteristic time scale of these clouds: 
R/AV. Clouds younger than this characteristic time didn't have 
time to be virialized and the virial theorem can not apply. This 
characteristic time scale for the SMC clouds is ~ 3 • 10 6 years. 
It is therefore possible that these clouds did not have time to 
be virialized. One possibility would be that, unlike Galactic gi- 
ant molecular clouds, in the SMC CO observations only traces 
clouds collapsing into stars on the timescale of a few million 
years (free-fall time). But this interpretation does not fit with the 
fact that two clouds, SMCB1#1 and Hodgel5, have no associ- 
ated Hn. 

In this study, v irial masses wer e deter mined using the sim- 
plified form from Mac Laren et al.l d!988l) . This form neglects 
among others the surface pressure and the magnetic field terms. 
The general virial equilibrium equation for a cloud in a steady 
state can be written as: 



= 2(7" - T ) + M + nV 



(5) 



(Mc Kee & Zweibelll992l) where T is the total kinetic energy 
(due to thermal and global motions), To represents the confine- 
ment by external thermal pressure, A4 is the magnetic energy 
and *W is the gravitational energy. The gravitational energy can 
be expressed as: 



3 GM 2 

<W = — a 

5 R 



(6) 



where G is the gravitationnal constant and a measure departure s 
fro m uniform density and sph ericity dBertoldi &McKee)ll99"2h . 
McK ee & ZweibeJ (119921) introduced the virial parameter 



a = a Jrnr so that the mass of a cloud is expressed as: 



M = 



5o- 2 R 
aG 



(7) 



If we take M\2mm as the true gas mass of the cloud, under- 
standing the high values of / = M\.2 mm jM v i r amounts to looking 
for low values of the virial parameters (with / = 10/9a), which 
is equivalent to solving: 



/ = 



10 [1-QPo/P)] 
9a [1 -{MIYW\)\ 



(8) 
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Table 5. Deduced characteristics for the SMC clouds. 



name 


./' 


MI\W\ 


\B\ 


n H 








liG 


10 J cm 3 


LIRS 49 


2.2 


0.5 


24 


0.8 


LIRS 36 


1.4 


0.2 


9 


0.3 


SMCB1#1 


2.0 


0.5 


12 


0.5 


SMCB1#2 


1.8 


0.4 


11 


0.3 


Hodge 15 


0.9 






0.4 


SMCB2 S 


3.1 


0.6 


20 


0.3 


SMCB2 N 


1.5 


0.2 


11 


0.4 



The values of f are then affected by density and geometry 
effects (a), external pressure (Pq) and magnetic fields (A1/| r U / |). 
As we saw in Galactic giant molecular clouds, / is smaller than 
unity and is interpreted as external pressure confinement. Since 
external pressure (high Pq) acts in lowering / and the effects of 
density and geometry are mild, we will only discuss the effects 
of magnetic field support in the SMC giant molecular clouds. 

5. 1 . Magnetic field support in clouds. 

From the observed characteristics of the molecular clouds, we 
can compute the magnetic-to-gravitational energy ratio needed 
to explain the / factors. Taking the case of a spherical uniform 
cloud with steady motions and a mean pressure exceeding the 
surface pressure (P » Pq), one obtains M/YW\ = 1 - ^. The 
obtained values are displayed in Table 0. 

Taking M = \{B 2 - B 2 )R 3 , one can express the strength of 
the magnetic field that is necessary to account for the / values: 



18 



" 5 \ 9f 



10 \ GM 2 



R 4 



(9) 



where Bo is the ambient medium value of the magnetic field 
strength. Taking Bq = 5/jG 0\ the magnetic field strength nec- 
essary to support the clouds is computed and displayed in Table 
[5] They range from ~ 10 to ~ 30 /uG, with a mean value of 
~ 15 /iG. Such stren gths are similar to the valu es reported for 
some Galactic clouds (iMvers & Goodmanl lT988) but the clumps 
surveyed are much smaller and confined b y the pressure, not 
by their self-gravity dBertoldi & M cKee 1992). Recently, a mag- 
netic field with a similar strength (~ 20yuG) h as been observed in 
the ou ter part of the Taurus molecular cloud (Wollebe n & Reichl 
12001 . 

So magnetic field support of the molecular clouds could ex- 
plain the mass difference observed in the SMC. But why would 
the magnetic field support be important for the SMC molecular 
clouds and not for the Galactic ones? Molecular cloud forma- 
tion could be different in the SMC. Giant molecular clouds in 
spiral galaxies are thought to form in spiral arms where the low 
shear and high gas density allows gravity to gather diffuse mat- 
ter over large scales. Matter may preferentially condense along 
field fines without comp ressing the magnetic field. We refer to 
Roso lowskv et al.l (|2003j) for a discussion of observational con- 
straints on proposed formation scenarios based on a survey of 
molecular clouds in M33. In an irregular galaxy like the SMC, 
massive stars have an important impact on the structure of the in- 
terstellar medium and molecular clouds may form out of shells 
blown by stellar winds and supernova e explosions. In the Spitzer 
infrared images dBolatto et al.l 120061) . the molecular clouds in 



5 consistent with lPohl il993h model; note that this choice does not 
really influence the result as long as the mean SMC magnetic field re- 
mains low, which we have no reason to doubt 



Fig. [TJdo appear as local peaks along a large shell of 300pc di- 
ameter. In this picture, the magnetic field is expected to be swept 
together with the matter. It is also possible that the lower metal- 
licity and therefore the lower dust content in the SMC creates 
lower dust opacity in the GMCs of the SMC than in our Galaxy. 
The ionization could therefore be higher in the SMC, supporting 
higher magnetic fields. Another possibility is that magnetic field 
support is important also in the Galactic molecular clouds and 
virial masses could have to be re-examined taking into account 
a magnetic field support on large scales. In this case, it would 
be interesting to see if the larger external pressure that will be 
required for the virial mass to match the other mass estimates 
(from y ray analysis for example) is compatible with observa- 
tions. Answering these questions are beyond the scope of this 
paper but will can be addressed in the future. In any case, mag- 
netic field observations (e.g through Zeeman effect) can be used 
to confirm the magnetic field strength in these regions. 

6. Conclusions 

We compared two different molecular gas tracers, the millime- 
ter dust emission and the virial theorem applied to CO obser- 
vations, at two different metallicities: our Galaxy and the SMC. 
SIMBA/SEST millimeter observations for the SMC and FIRAS 
observations for the Galaxy have a similar linear resolution, en- 
abling a uniform study on similar giant molecular clouds in these 
two environments. 

We find a clear difference between the two samples: in the 
Galaxy, the virial mass is systematically larger than the mass 
deduced from millimeter emission, while in the SMC sample, 
the virial mass is lower than the mass deduced from millime- 
ter emission. The main uncertainty in the mass determination 
comes from the dust emissivity per hydrogen atom in the SMC. 
We have used the millimeter emissivity determined in Galactic 
molecular gas scaled by 1/6 to account for the lower SMC metal- 
licity. This scaling factor has been calculated asssuming similar 
depletions of heavy elements in the SMC and in the Galaxy. We 
consider this assumption to be conservative because diffuse ISM 
depletions have been inferred to be smaller in the SMC than in 
the Galaxy from UV spectroscopy and dust observations. The 
millimeter cloud masses would have been larger for smaller de- 
pletions (i.e. smaller dust-to-gas ratio), enhancing the discrep- 
ancy. 

For Galactic clouds the dust and CO masses are consistent 
with the gamma ray detemination of the Galactic Xqo factor. 
The virial masses are higher than the cloud masses. The under- 
estimation of SMC molecular cloud masses by virial methods is 
thus an unexpected result. Inspecting the assumptions made in 
the virial mass computations, we find that magnetic field sup- 
port in the clouds could explain the difference between the two 
mass estimates. In this scenario, SMC giant molecular clouds are 
not fully supported by turbulent motions. Additional support is 
given by a magnetic field of ~ 15/jG Further studies and mag- 
netic field observations are needed to confirm this interpretation. 

An alternative explanation is that the virial mass is under- 
estimated because the CO line width does not measure the full 
amplitude of the turbulent motions within the cloud. In the SMC 
GMCs, the CO emission is expected to come from dense clumps 
embedded in lower density gas where CO is photo-dissociated. 
The bulk of the cloud mass is in the lower density gas, but the CO 
line width only measures the velocity dispersion of the clumps. 
Can the lower density gas turbulence be higher than the clump to 
clump velocity dispersion? This study emphasizes the power of 
millimeter emission as a dense matter tracer. The large amount 
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of millimeter data coming will enable interesting follow up stud- 
ies with better resolution and sensitivity. 
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Appendix A: Emissivity of dust associated with 
molecular gas 

Dust emissivity per hydrogen atom is often taken in the dif- 
fuse medium where it is easy to deduce directly from ob- 
servations. Howev er, dust properties cou l d vary from diffuse 
to dense region s (Cambres y et alj 1200 it IStepnik et alj [2003; 
ICambresv etaT] l2005h . We therefore chose to recompute the 
emissivity of dust is computed in the molecular ring in order 
to represent the properties of dust associated with molecular hy- 
drogen. 

In the molecular ring, the 1.2mm dust continuum and CO 
emissions longitude profiles along the Galactic plane are well 
correlated and the ratio of their intensities is constant: R = 
(0.89 + 0.01) ■ W- l MJy(K.km.s- l .sry l (Boulanger et al., in 
prep.). This ratio is related to the well known CO-to-H2 con- 
version factor (X = N(H2)I Ico) through the following equation: 



R = 2e H (A)B A (T dusr )X, 



(A.l) 



where en is the dust emissivity per hydrogen atom in the ob- 
served region. Although Td ust may vary from one cloud to the 
other, these variations are likely to be moderate and can therefore 
be neglected at millimeter wavelenght while T > 10K. In this 
study, we consider the temperature range measured by FIRAS as 
pertinent for molecular clouds analysis. The CO-to-H2 factor is 
empirically determined in different environments and is largely 
discussed in the litterature. 

In our Galaxy and in particular i n the solar neigh bourhood, 
studies of y rays and gas tracers (Stron g~et alj 1 19881) enable 
a determination of this X factor. For our Galaxy as a whole, 
the analysis of recent observations gives X factors between 
1.56 and 1.9 X \tf a molcm" 2 .{K.km.s^Y l dStrong & Mattoxl 
119961: iHunteret all 1 19971: iGrenieret al.ll2005l) . We use the re- 
cent estimate obtained for the local molecular clouds through 
y rays analysis (Grenier et al., priv. comm.): X = 1.8 + 0.3 • 
lO^mol.cm^iK.km.s- 1 )- 1 . 

Using this local Galactic X factor and using the ratio of mil- 
limeter and CO intensities determined in the molecular ring, we 
compute the emissivity of dust associated with molecular gas in 
our Galaxy: 



e™ f (1.2mm, H 2 ) = (1.20 ± 0.20) x lQ- Ib ar L .cm 



26„.-l 



(A.2) 



This value is used in eq.|2]to compute molecular cloud masses. 
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The dust emissivity in the diffuse medium is measured 
through the correlation at high Galactic latitude betwwen 
dust emission and H i emission at low column densities 
(Boulange r et al.|[l996l) . For a dust temperature of 17. 2K, this 
correlation gives a dust emissivity per hydrogen atom of: 

e% f (12mm,HI) = (6.7 + 1.2) • 10~ 21 at' 1 . cm 2 . (A3) 

Dust emissivity in the molecular medium is therefore twice 
larger than in the diffuse medium. Such an increase in emis- 
sivity from the diffuse medium to molecular clouds was ob- 
served in our Galaxy, but on much smalle r and denser regions 
dCambresv et al.ll200ltlStepnik et a l. 2003). This increase is ob- 
served to be correlated with a drop small grain abundance. It has 
been thus interpreted as the coagulation of small grains between 
themself or on large grains in dense regions of molecular clouds, 
giving composite porous grains with a larger emissivity. 



Comparing this value to the dust-to-hydrogen mass ratio for £ 
Oph dus(J the dust-to-gas mass ratio difference between the 
SMC and the Solar Neighbourhood is found to be: 



Appendix B: The dust to gas ratio in the SMC 



The gas-to-dust ratio in the SMC has been deter mined by us- 
ing extinction curves to derive N(H i)/A(V) values dGordon et al.l 
2003: lBouchet et alJll985MTumlinson et al.ll2002l) . However, the 
derived values correspond to specific line of sights (where ex- 
tinction can be measured), they rely on the present extinction law 
understanding, and there is a significant scatter f rom one region 
to the other and from one study to the other R iBot et all 12004) 
deduced a gas-to-dust ratio 30 times the solar value by compar- 
ing the dust and H i emission of the SMC. However, this value 
refers to the diffuse medium and could reflect dust grain destruc- 
tion processes that do not apply so intensively in the denser en- 
vironnements of the SMC bar where grains are more shielded. A 
generally used value for the SMC dust-to-gas ra tio is one tent h 
of the solar one, using the metallicity difference (lDufour| [l984). 
However, it relates to gas phase abundances while we are inter- 
ested in dust masses and most im portantly since 1984 the solar 
abundances have been updated (Sofia & Mever 2001bja|). 

We therefore choose to recompute the SMC dust abundances 
to evaluate the dust-to-gas mass ratio difference between the 
SMC and the Solar Neighbourhood. For that, we assume that 
the elemental depletions are the same in the SMC as in the Solar 
Neighbourhood. This is supported by studies of Sk 108 in the 
main bar of the SMC, where the depletion patte r ns are similar to 
those o f Galactic halo clouds lWeltv et al.l 0997); Mall ouris et all 
d200ll) . The depletions in the SMC could be less than the solar 
values dRolleston et al.l 12003). so that our estimation could be 
an upper limit. Comparing the ( Oph dust and solar abundances 
given by I Sofia & Mever! (2001b a), we inferred depletions r\ of 
0.66, 0.1 1, 0.46, 0.92, and 0.95 for C, N, O, Mg, and Si respec- 
tively. 

The main sequence star AV 304 shows no sign of process- 
ing in its photosphere and could therefore repres ent the present- 
day to tal (gas and dust) compositio n of the SMC (Roll eston et al.l 
2003). Using lTrundle et a l. (2004) corrected abundances for this 
star with the depletion factors computed above, the dust-to- 
hydrogen mass ratio can be computed: 



= y [M/ H] t0 , ■ T]m — , 

m H tjf m H 



(B.l) 



and is used as an estimate of the dust-to-gas mass ratio. We ob- 
tained a dust-to-hydrogen mass ratio of 1 .5 1 ■ 10~ 3 for the SMC. 



6 The dust to has ratio found this way to range between ~ 1/6 and 
-1/30 the Galactic value 



7 nidi st/>nH((Oph) = 8.88 • 10 3 as computed from ISofia & Meveil 
fcOOlbl lal) abundance values. 



